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Abstract. 

The reconstruction scheme is developed for modified f(R) gravity with realistic 
matter (dark matter, baryons, radiation). Two versions of such theory are constructed: 
the first one describes the sequence of radiation and matter domination, decceleration- 
acceleration transition and acceleration era and the second one is reconstructed from 
exact ACDM cosmology. The asymptotic behaviour of first model at late times 
coincides with the theory containing positive and negative powers of curvature while 
second model approaches to General Relativity without singularity at zero curvature. 
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1. Introduction 

The modified f(R) gravity (for a review, see [Tj) suggests very interesting gravitational 
alternative for dark energy where the current cosmic speed-up is explained by the 
presence of some sub-dominant terms (like 1/R [21 [H] or InR [I]) which may be caused by 
string/M-theory [SJ. These terms may become essential at small curvature. There are 
simple models of modified f{R) gravity like the one with negative and positive powers 
of curvature [0] which are consistent with late-time astrophysical constraints IB] and 
solar system tests jU]. The accelerating cosmology in such theories has been studied in 
rets.© EH]. 

It became clear recently that modified f{R) gravity may also describe the sequence 
of matter dominated, transition from decceleration to acceleration and acceleration eras 
[TT| . It has been developed the reconstruction scheme which permits to reconstruct 
modified gravity for any given FRW cosmology JJJ. In the present communication we 
show that such reconstruction scheme may be applied also in the realistic situation when 
usual matter is added to modified f(R) gravity. Specifically, we find two implicitly 
given versions of f(R) gravity with matter which may serve as an alternative for 
ACDM cosmology. The first model describes the sequence of radiation dominated, 
matter dominated and acceleration eras while second model matches exactly with 
ACDM cosmology. Moreover, in the acceleration epoch the asymptotic behaviour of 
such theories may be defined: in first case the model [H] containing negative power 
of curvature is recovered while second model asymptotically approaches to standard 
General Relativity. 

2. General formulation of the reconstruction scheme 

In the present section we review the reconstruction scheme for modified gravity with 
f(R) action ^JJ, where it has been shown how any given cosmology may define the 
implicit form of the function /. The starting action of modified gravity is: 

S = J d 4 Xy^f(R) . (1) 

The above action is equivalently rewritten as (see, for example, |13j ) 

S = J d^^ {P(<J>)R + Q{(j>) + Matter} ■ (2) 

Here P and Q are proper functions of the scalar field <fi and £ m attcr is the matter 
Lagrangian density. Since the scalar field does not have a kinetic term, it may be 
regarded as an auxiliary field. In fact, by the variation of <ft, it follows = P'(<f))R+Q'((j)), 
which may be solved with respect to (ft as (ft = <p(R). By substituting the obtained 
expression of <f)(R) into (J2J), one comes back to /(-R)-gravity: 

S = J >aV^{/0R) + /: mattcr } , f(R) = P(<f ) (R))R + Q(<P(R)) . (3) 
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By the variation of the action (J2J) with respect to the metric g^, we obtain the equations 
corresponding to standard spatially-flat FRW universe 

= - w 2 P(4>) - Q(4>) - m^P- + p , (4) 

= (AH + m + QM + 2^™ + 4ff « +p . (5) 

Simple algebra leads to the following equation [TT] 

d 2 P((b(t)) dP(6(t)) ■ , s 

= 2 ^>l-2H K j} +AHP )+p+p. 6 

at^ d<p 

As one can redefine the scalar field <fi freely, we may choose <fi = t. It is assumed that 
p and p are the sum from the contribution of the matters with a constant equation of 
state parameters u>j. Especially, when it is assumed a combination of the radiation and 
dust, one gets the standard expression 

-4 | -3 PrO -4 fry\ 

p = p r0 a + p do a , p = —a , (7) 

with constant p r Q and pdo- If the scale factor a is given by a proper function g(t) as 
a = a e 9 ^ with a constant a , Eq.(jHJ reduces to the second rank differential equation 
(see also [T^]): 

+ £(1 + «*) Pi0 ao 3(1+ ^e- 3 ^) ff W . (8) 

i 

In principle, by solving (jHJ) the form of P{4>) may be found. Using (J3J) (or equivalently 
(JSJ), the form of Q(0) follows as 

W) = "6 (^(0)) 2 P(0) - 6^)^ +EAoao 3(1+Wi) e- 3(1+ ^ )s( ^ -(9) 

Hence, in principle, any given cosmology expressed as a = aoe 9 ^- 1 can be realized as the 
solution of some specific (reconstructed) f(R) -gravity. Moreover, this reconstruction 
scheme as is shown in ref. [13] may be generalized for other types of modified gravity. 



3. Models of f(R) gravity with transition of matter dominated phase to the 
acceleration phase 

Let us consider realistic examples where the total action contains also usual matter. 
The starting form of g(<p) is 

s(0) = M0)ln(£) , (10) 

with a constant 0o- It is assumed that h(4>) is a slowly changing function of 0. Due to 
(j) — t, it follows H ~ h(t)/t, and the effective EoS (equation of state) parameter w e s is 
given by 

2H 2 
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Therefore w e g changes slowly with time. 

As h(<p) is assumed to be a a slowly changing function of <fi — t, one can use 
adiabatic approximation and neglect the derivatives of h(<p) like {h'(<f)) ~ h"(<f)) ~ 0). 
Then Eq.(JHJ) has the following form: 

o = M_l«_^ )+E (12) 

d(p z <p dt <p z i 

The solution for P(4>) is found to be ^1] 

P(0) = + v _^- W + J2 p l ((j))(j)~ 3il+m)hW+2 . (13) 

i 

Here p± are arbitrary constants and 

M0)-i±^/M0) 2 + 6M0) + i 

n±{<t>) = —2 . 

Pi(<i>) = -{(l + w^oa^^Ul^^} 

x {6(1 + Wi )(A + Zw t )h((j)) 2 - 2 (13 + 9wi) h{4>) + 4p . (14) 
Especially for the radiation and dust, one has 

PradK</>) = - 



3a£ (4O/i(0) 2 - 32/i(0) + 4) ' 

Pdust[(p) - (24/i(0) 2 - 26/i(0) + 4) ' 1 J 

The form of Q(</>) is found to be 

Q{4>) = - 6h{<l>)p + (h(<j>) + n + (0)) n +^- 2 
- 6/i(0)p_ (/i(</>) + n_(0)) n " w " 2 
+ E {-6/^(0) (-(2 + 3 Wi )/i(0) + 2) Pi {4>) 

i 

+fto«o 3(1+ " l Vo (1+U,l) " W }r 3(1+ ^ )M</,) ■ (16) 

Eq.fllUJ) tells that 

6 (-h(t) + 2h(tf 

T 2 

Let assume lim^oM^) = ^ an d lim^oo /i(</>) = ft,/. Then if < hi < 1, the early 
universe is in decceleration phase and if /i/ > 1, the late universe is in acceleration phase. 
We may consider the case h(4>) ~ h m is almost constant when ~ t m (0 <C t m <C +oo). 
If Tii , /i/ > 1 and < h m < 1, the early universe is also accelerating, which could 
describe the inflation. After that the universe becomes deccelerating, which corresponds 
to matter-dominated phase with h(<f)) ~ 2/3 there. Furthermore, after that, the universe 
enters the acceleration epoch. Hence, the unification of the inflation, matter domination 
and late-time acceleration is possible in the theory under consideration. 

As an extension of the above model JI] , we consider the inclusion of the radiation, 
baryons, and dark matter: 

m = ~TTo¥~ ' (18) 



T 2 • ( 17 ) 
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with constants hi, hf, and q. When — > 0, h(<f)) — > hi and when — > oo, /i(0) — > /ij. If 
g is small enough, h(<f)) can be a slowly varying function of 0. By using the expression 
of (H3), we finddU 

2 = $ (i?) , $ ± (i?) , 

$ = a V3 + a V3 $ = «V 3 e 2 -/ 3 + a^V 2 ™/ 3 , 



/So 



2 

2 (2i? + 6h f q - 12/iJg 



27g 3 i? 3 

2i? + 6h f q - 12/iJg) (i? + 6h iq + 6% - Ahih f q) 

3q~R 

+ 6hi - 12/i 2 , 

_ (2R+6h f q- 12%) 2 R + Qh iq + 6h f q~Ahih f q 
A = 3?^ • (19) 

There are three branches $o and $± in (j!9jh Eqs. ()17j) and (|T8|) show that when the 
curvature is small (0 = t is large), we find i? ~ 6 ft,/ + 2/i 2 ) /0 2 and when the 
curvature is large (0 = t is small), ~ 6 (—/?-« + 2/i 2 ) /0 2 . This asymptotic behaviour 
indicates that we should choose $o i n Q19|) . Then explicit form of f(R) could be given 
by using the expressions of P{4>) (O and Q(4>) (|TB)l as 



f(R) = p (yJMR)) R + Q [yJMR)) ■ (20) 

One may check the asymptotic behavior of f(R) (|2(Jj) (for some parameters choice) 
in the acceleration era coincides with the theory proposed in 6J. We now consider the 
case that, besides dust, which coul be dark matter and baryons, there is a radiation. In 
this case, P{<j>) is given by 

P(0) = p + n +W + (p dark (0) + p baryon (0)) 0" 3 ^)+ 2 + p rad (0K 4M0)+2 (21) 

Here p rad (0) is given by flE} and p dar k(0) + Pbar yO n(0) are 



IPdarkO + PbaryonO ) <Pq 

~a~l (24/i(0) 2 - 26/i(0) + 4) 



Pdark(0) +Pbaryon(0) = 3 /„,, / ,x 2 ngTTXN I /A ' ( 22 ) 



We now find n + > -3/i(0) + 2 > -4/i(0) + 2 > 0, in <f5Tjl - Here n + = 
(/i(0) - 1 ± ^//i(0) 2 + 6/i(0) + l) /2 is defined in ([14)1 . Then when is large, the first 
term in (|2*Tj) dominates and when is small, the last term dominates. When is 
large, curvature is small and 2 ~ 6 (—hf + 2hf) /R and h(<f) — > h(oo) = hf. Hence, 
f(R) ~ jr(^- H \/^P)/ 4 . Especially when /i » 1, it becomes /(#) ~ R~ h fl 2 . 
Therefore there appears the negative power of R predicted by the presence of matter- 
dominated stage. As H ~ hf/t, if hf > 1, the universe is in acceleration phase. 

On the other hand, when curvature is large, we find 2 ~ 6 (—hi + 2hi) / R and 
h(<f>) — > h(0) = hi. If the universe era corresponds to radiation dominated phase 
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{hi = 1/2), P(4>) becomes a constant and therefore f(R) ~ R, which reproduces the 
Einstein gravity. 

Thus, in the above model, the radiation/matter dominated phase evolves into 
acceleration phase and f(R) behaves as f{R) ~ R initially while f(R) ~ 
jj-(h(<j>)-5+y/h 2 f +6hf+i)/4 ^ j a ^ e ^ me p or foe matter dominated phase, we have h = 2/3. 
Since hi = 1/2 < 2/3 < hf and h{4> — t) is a slowly increasing function of = t, there 
should always be a matter dominated phase. Therefore in the model ()18j) with hi = 1/2, 
the universe is first in radiation dominated phase. Subsequently, the universe evolves to 
the matter dominated phase, and finally to accelerated phase which is consistent with 
ACDM cosmology. 

Thus, we presented the example of implicitly given f(R) gravity which describes 
the radiation dominated era, the matter dominated stage, transition from decceleration 
to acceleration and acceleration epoch (where it may include the negative powers of R). 
This model seems to be quite reasonable alternative for the standard ACDM cosmology. 



4. Model reproducing ACDM-type cosmology 

Let us investigate if ACDM-type cosmology could be reconstructed exactly by f(R)- 
gravity in the present formulation when we include dust, which could be a sum of the 
baryon and dark matter, and radiation. 

In the Einstein gravity, when there is a matter with the EOS parameter w and 
cosmo logical constant, the FRW equation has the following form: 

^H 2 = Poa~ 3(1+w) + 4t2 ■ ( 23 ) 

K z K z l z 

Here I is the length parameter. The solution of (|2*3*|) is given by 

a = aoe**) , 9(t) = In (a sinh (f - to)) ) ■ (24) 

Here t is a constant of the integration and a 2 = (l/3)K 2 l 2 p a ~ 3( ' 1+w \ Let us show 
how is possible to reconstruct /(-R)-gravity reproducing (J2H). When matter is included, 
Eq.fjHJ) has the following form: 

n n d 2 PU) 2 , /3(1 +w) , , ,\dPU) 
= 2-^- y coth 



6(1 + w) . , 2 I 3(1 + w) A „. A 

1 p } sum" 2 f V % > (0 - t ) j P(<P) 



4 _ 4 / . (Z{l + w) 
H — Pro^o ot smh 1 



3 r u V V 21 



-8/3(1+™) 



-2/(l+w) 



+ Pdo a 3 L sinh (0-t o )J ) • (25) 

Since this equation is the linear inhomogeneous equation, the general solution is given 
by the sum of the special solution and the general solution which corresponds to the 
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homogeneous equation. For the case without matter, by changing the variable from 
to z as follows, 

^-sinh- 2 (^±^(t-t )) , (26) 

Eq. (|25j) without matter can be rewritten in the form of Gauss's hypergeometric 
differential equation: 



dP 



7 - (a + $ + l) z] af3P 



dz 

7 = 4 + Qnl v 5 + ^ + 1 = 6 + Qnl v ^ = - on! v ^ 
3(1 + w) 3(1 + w) 3(1 + w) 

whose solution is given by Gauss's hypergeometric function: 

P = P i<(a,^,7, z) = P — — — 2^ — v r • I 28 ) 

Here T is the T-function. There is one more linearly independent solution like 
(1 — z)^~ a_|8 F(7 — a,7 — (3, 7; z) but we drop it here, for simplicity. Using (j^J), one 
finds the form of Q(4>)'- 

Q= - 6{1 ~/ )Po F(aJ,r, z ) 

From (J26j) . it follows z —>■ when t = </> — > +oo. Then in the limit, one arrives at 
P(0)P + Q{(t>) -> P P - 6P /Z 2 . Identifying P = 1/2k 2 and A = 6// 2 , the Einstein 
theory with cosmological constant A can be reproduced. The action is not singular 
even in the limit of t — > oo. Therefore even without cosmological constant nor cold 
dark matter, the cosmology of ACDM model could be reproduced by /(i?) -gravity (for 
a different treatment, see [TH]). 

We now investigate the special solution of (|25|). By changing the variable as in (|26|). 
the inhomogeneous differential equation looks as: 



d 2 P 



7 - (a + p + l) z] ^ - a/3P 
v 7 J dz 



l n-2(1+3w)/3(1+™) , e / n-t^ 
A/ 2 7 2 

77 " 27(1 + W ) Pr ° G ° " ' C " 3(1 + w) Pd ° " • (dUj 
It is not trivial to find the solution of ([30)1 . Let us consider the case that w = and 
z — ► — oo, that is, £ — ► t - I n the limit, Eq.(j3Ul) reduces to 

= -z 2 ^ + 72^ - afiP + r] (-z)- 2/3 , (31) 
dz z dz 

whose special solution is given by 

P = P (-,)- 2 / 3 , P = = -%■ (32) 
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In principle, there could be found other special solution of Eq. (fHUj). This proves that 
even in the presence of matter, the standard ACDM cosmology could be reproduced by 
/(-R)-gravity exactly 

Thus, we presented two versions of modified f(R) gravity with matter. The first 
version describes the sequence of radiation dominated, matter dominated, transition 
from decceleration to acceleration and acceleration eras (compare with scalar-tensor 
gravity with the same emerging cosmology JE]). I n the aceleration era the action 
may asymptotically approach to the action with negative and positive powers of R 
proposed in 6j. Moreover, for some choice of parameters it reproduces the ACDM 
cosmology at late times. The second model may reproduce ACDM cosmology exactly. 
Using the number of free parameters of the models one can expect that they may be in 
good correspondence with observatonal data as they are with three years WMAP data. 
Nevertheless, the precise fitting of the proposed f(R) gravity against existing/coming 
observational data should be done. 
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